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1. Introduction 



Polynomials are fascinating because they have many facets to their personalities. By 
definition, a polynomial / G C[x] is an expression of the form 

(1.1) f(x) = a + a\x + ■ • • + a n x n , 

where the cij are complex numbers. By the fundamental theorem of algebra, / also has a 
representation as a product 

(1.2) f(x) = a n (x - x x ) ■ ■ ■ (x - x n ), 

where Xj are the roots of /. The fact that there are two ways of looking at polynomials 
provides possibilities that are hidden by their apparent simplicity. 

In this paper we look at two different polynomials which have equally spaced zeros, and 
we study the zeros of their derivatives. We will see that Euler's famous result 

• ' 00 i 2 

S : (!'3) C(2):=E^ = T 

■rr: ■ n z 6 

OO . n=l 

is lurking in the background, as are some generalizations of Euler's result. Our proof that 
C(2) = 7t 2 /6 appears to be new, and it doesn't use anything more than is typically covered 
^ ■ in a second semester calculus course. 

2. A Polynomial With Zeros at the Integers 



The first polynomial we will consider is the degree N +1 polynomial with simple zeros at 

0. 1.2, ...,N. That is, 

N 

(2.1) p N (x) := x(x - 1) ■ ■ ■ (x - N) = Y[(x - n). 

n=0 

We will write p(x) instead of Pn{x) to simplify the notation. The question we will obtain a 
partial answer to is: Where are the zeros of the derivative p'(x)7 

By Rolle's theorem, we know that p' has a zero in each of the intervals (0, 1), (1, 2), . . . ,(iV— 

1, N). Those zeros account for iV of the zeros of p', which is all of the zeros because p' has 
degree N. 

We are left with the question: where in those intervals (0, 1),. . . ,(N — 1, N) are the zeros of 
p'l Suppose more generally that / is a polynomial with only real zeros. There is a common 
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misconception, reinforced by looking at small degree examples and badly drawn graphs, that 
the zeros of /' should be located close to the midpoint between neighboring zeros of /. A 
more accurate guide is the following: 

Repulsion Principle. Zeros of f try to move away from the zeros of f . 

What the principle says is: when /' decides where to put a zero, it looks around to see 
what zeros of / are nearby If it sees a bunch of zeros of / to the right, it will try to shift 
the zero of /' to the left, and vice-versa. This is illustrated in Figure [I] which shows part of 
the graph of a high-degree polynomial. 




Figure 1. The circles are the midpoints of neighboring zeros of the polynomial /, 
and the dots are the zeros of /'. The zeros of / push the zeros of /' away from 
the midpoints. 



Let's focus on the interval (0, 1) for the polynomial p(x) given above. Since there are 
many zeros to the right, and only one zero to the left, the Repulsion Principle says that p' 
will shift the zero to the left, so the zero will be closer to than to 1. The question is: how 
close? We answer that question in the next section. 

The Repulsion Principle, not by that name, is discussed extensively in our paper [TJ. This 
principle explains, for example, the observation in [I] that the smallest zero of /' is closer to 
the smallest zero of / than to the next-smallest zero. The specific case of the polynomial pn 
is discussed in [5]. 

3. Finding zeros of p' 

Any attempt to find a simple formula for the location of the first zero of p' will fail for 
two reasons. First, there is no simple formula for the derivative p' . We discuss this difficulty 
in the next section. Second, there is no closed form expression for the zeros of an arbitrary 
high degree polynomial, so even if we had a nice formula for p'(x) it might not do us any 
good. 

3.1. A formula for p'(x). To see the difficulty of finding a formula for p'(x), first consider 
the product representation (12. ip . To differentiate a product of N + 1 things, you have to use 
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the product rule N times. The result will be a sum of N + 1 terms, where each term is an 
iV-fold product. The entire expression is quite unwieldy. 

Since differentiating a sum is easier than differentiating a product, we can try multiplying 
out (12. f I) to write p(x) as a sum: 



(3-1) p(x) = 



N+l 

d n X 
n=0 



A few of the coefficients are easy to find. For example, ao = because p(0) = 0, which 
we see by substituting x = in the product. It is not too hard to see that a\ = (—1) N N\ 
and ajv+i = 1. A couple more coefficients can be described as simple formulas, but there is 
no known closed formula for all of the coefficients. These coefficients come up in a variety 
of situations and are know as the signed Stirling numbers of the first kind. We have a n = 
s(N + 1, n) in the usual notation for the Stirling numbers. Since, there is no simple formula 
for s(N + l,n), differentiating (13. f I) does not give a useful expression. For more information 
about the Stirling numbers, see Chapter 6 of [2]. 

Since the product form of p is difficult to differentiate, and the sum form is easy to 
differentiate but has unusable coefficients, what can we do? The answer is: take logarithms! 
Since the log of a product is a sum, we have 

log(p(x)) = log I - n) J 

\n=0 / 

N 

(3.2) = ^log(x-n). 



n=0 



Note: here "log" means natural log, the function that calculator manufacturers and non- 
mathematicians call "In". 

Differentiating the above expression we obtain "the logarithmic derivative" of p: 

(3.3) 5^7 = £_*_. 

p[x) x — n 

Fy > n=0 

This is useful because the solutions to p'{x) = are the same as the solutions to 

(3.4) ^ = 0, 

p(x) 

except possibly at the points where p(x) = 0. However, we are trying to understand the 
solution to p\x) = which is between and 1, so we have the following equivalent problem: 
find the solution to 

N 1 

^x-n ° 

n=0 



which satisfies < x < 1. We will not find an exact expression for that solution, but we will 
find a good approximation. 
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3.2. First approximation to the solution. In this section we work informally and make 
things more precise in the following sections 

Let a denote the solution to p'(x) = with < a < 1. As described in the previous 
section, we must solve 

(3.6) = - + - + ••• + 



a a — 1 a — N 

Since solving (13. 6p is the same as finding a root of an iVth degree polynomial, we should not 
expect to find an exact formula for a, except maybe for small N. Instead, we will look for 
an approximation. 

Recall that the Repulsion Principle told us that a is close to zero, so in ( 13. 61) the term 
1/a is much larger than the rest. So we will rearrange (13.61) to put 1/a by itself: 

(37) i = + — !— + •••+ 1 



a 1 — a 2 — a N — a 

Since a is close to 0, the right side of f ]3. Tl) will be close to the Nth harmonic number: 

(3-8) + 1 + + 

We give more information about the harmonic numbers in the next section, but for now just 
recall that is approximately logiV. Thus, 

(3.9) — « log N. 

a 

so, a ~ 1/logiV, and we see that, for large N, a is actually quite close to 0. In particular, 
a is not close to the midpoint between the first two zeros. 

The above calculation is informal and we do not intend to convey a precise meaning with 
the "approximately" symbol "~" . We will do the calculation rigorously in the next section. 

4. Better approximations to the first zero of p' 

We make precise the argument in the previous section, and find that it touches on many 
combinatorial and number-theoretic objects that have been widely studied. In our approxi- 
mation to the first root of p'{x) we will encounter Euler's constant 7 and the Riemann zeta 
function 



k=l 

Then in Section [5] we will differentiate another polynomial in order to evaluate C(2), and in 
Section [6] we evaluate the zeta function at the even integers. That calculation will involve 
the Bernoulli numbers, the Stirling numbers, and the Eulerian numbers. 

To give a precise estimate for < a < 1, the first root of p', we need two facts. First, the 
geometric series 

(4.2) = 1 + x + x 2 + ■ ■ ■ + x M + 0(x M+1 ) } 

1 — x 

which can be manipulated to give 
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as x — > 0. Second, 

N 1 f l\ 

(4.4) Hn :=J2- = log(^V) + 7 + O ( ^ ) • 

n=l 71 ^ ' 

Here 7 is Euler's constant, 7 ~ 0.577. Equation (14.41) is the definition of Euler's constant, but 
it appears in many other places, such as the Laurent series for the Riemann-zeta function: 

(4.5) C(s ) = _l_ + 7 + 0( s -l). 

There is a whole book on this interesting number [3]. 
We have 

Theorem 4.1. I/O < a < 1 is a root of p' N (x) then 

(4.6) - = log(iV) + 0(l), 

a 

so 

<47) a ^ + °yw 

Proof. The only fact we will use about a is that < a < 1, so in particular a = 0(1). Start 
with equation (13. 7p . and use equation (14. 3 p with M = 0, and (14. 4D : 

1 N 1 

1 = E— 

i=i J 



^ /1 



= log(iV) + 7 + ^ ( jy ) + since 1/n 2 converges 

= log(JV) + 0(l). 

That gives the first equation in the theorem. To obtain the second, take the reciprocal of 
both sides and use (14. 3[) with j = log N and x = 0(1). □ 

We started with the assumption < a < 1 and found an approximation for a with an 
error term 0(1/ log 2 (N)). By feeding this into an almost identical calculation, we can obtain 
an even better approximation to a. 

Theorem 4.2. If < a < 1 is a root ofp' N (x) then 
so 

(4-9) « = ^-^ + ?4i + 



log(iV) log 2 (iV) log 3 (iV) Vlog 4 (iV) 
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Euler proved that 

00 1 2 

(4-10) C(2) := X) =5 = T- 

n=l 

In the proof below we will only use the fact that 

n=l n=l n=N+l x ' 

which follows from 

1 C 00 1 °° 1 f 00 1 1 

(4.12) = / — dx < V — < / — = — . 

K J N + l J N+1 x 2 n ^ +1 ™ 2 Jn N 

Proof of Theorem \4.S\ Again we start with equation (13. 7p and use equation (14. 3p . this time 
with M = 1, and (@3D- We have 

iV 



3=1 J 
* / 1 



iV /, / 2 

j=1 \3 3 \3 3 

n n N /^\ 

3=1 J 3=1 J 3=1 VJ 7 

log(iV) + 7 + O ( + a (((2) + O ( i ) ) + Oi < r ) 



The last step used the fact that converges. On the final line above, replace a by 

1/ logiV + 0(1/ log 2 N) and 0(a 2 ) by 0(1/ log 2 N), both of which follow from TheoremEJ 
Since 1/N = 0(1/ log N) this establishes the first formula in Theorem 14.21 The second 
formula follows in exactly the same way as the second formula in Theorem 14. 11 □ 

It is clear that we could use the result in Theorem I4.2[ feed it back into equation (13. Tj) , 
and obtain an even more precise formula for a. Repeating this process would establish that 



(413) a -io^ + k^ + --' + io^ +0 n^ 



where the Cj are polynomial expressions in 7 and C(2), ■ • • , C(j — !)• F° r the first few Cj we 
have 

c 2 = - 7 

c 3 =7 2 - C(2) 

c 4 =C(3)-7C(2)+C(2)-7- 

In the next section we use a different polynomial, but which also has equally spaced zeros, 
to prove Euler's result £(2) = 7r 2 /6. 
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5. The polynomials z n — 1, and ((2) 

The previous section we considered polynomials with equally spaced zeros on the real line. 
In this section we consider polynomials with equally spaced zeros on the unit circle. That 
is, q(z) = qN{z) '■= z N — 1. If we set e(z) := e 2mz then the zeros of q(z) are at the iVth roots 
of unity: 1, e(l/N), e(2/N),. . . ,e((N - 1)/N). Thus, 



(5.1) z N -l= J] (z-e(j/N)). 

j mod N 

Here j mod N refers to a set of representatives of the integers modulo N. This is legitimate 
because e(z) has period 1. Later we will use the specific choice —N/2 < j < N/2. 

We now use the polynomial q(z) to evaluate C(2)- As in the previous section, we begin 
with the logarithmic derivative of q. Since we are trying to prove an identity, in this case 
we will make use of both the product and the sum representation of the polynomial. Taking 
the second derivative of the logarithm of (15.11) we obtain 

—N(N - l)z N - 2 - Nz 2N - 2 _ x ^ 1 

(z N - I) 2 ~~ .^i„ (z-e(j/N)) 2 ' 

y ' j mod N y w ' " 

Various formulas can be obtained by plugging in specific quantities for z in the above equa- 
tion. It is tempting to plug in z = 1, which is possible if one first moves the j = 
mod N term to the left side and simplifies. Equivalently, start over with the polynomial 
(z N — l)/(z — 1). Instead, we will plug in z = e(l/2N), a 2iVth root of unity, which is 
motivated by noting that the left side will become fairly simple. 

Multiplying both sides of that equation by z 2 and plugging in z = e(l/2N) we have 

Lemma 5.1. We have 

N 2 -2N x ^ 1 



-e(l/A0 Y, 



, {X AV^)-e{j/N)f 

(5.3) =- y x. 

(1 - e((2j - l)/2iV)) 2 

The remainder of the calculation closely parallels the treatment in the previous section: 
we want to expand the denominator of the right side of (15.31) as a series. Using the Taylor 
series for e x one can check that for < \x\ < 1 we have 

(5.4) - = — + — ^— + — - -inx + 0(x 2 ). 

K 1 (l_ e (x)) 2 4vr 2 a; 2 2tix 12 6 V ; 

The restriction on x is due to the fact that e(x) = 1 for integral x. We will only need the 
first term in that expansion. 

Since (15.41) is only valid for < \x\ < 1 we are forced to make the specific choice —N/2 < 
j < N/2 for the representatives of the integers modulo N. Making this choice and then 
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applying the Taylor expansion (15.41) we obtain 



N 2 - 2N _ v ^ 1 

~ ' -«k<m {l - e ^- l) l 2N))2 

^ V4vr 2 (2j-1) 2 V 2? + 1 

-N/2<j<N/2 V V J ' 

N 2 1 



, 2 E j^^ + o^n). 

-N/2<j<N/2 v J ' 

Dividing both sides by iV 2 and rearranging, we have 

(5 - 5) ^ (2j^ = T + hv-J- 

-N/2<j<N/2 v J y v 7 

Letting iV — > oo and replacing j — > 1 — j for j < 0, we obtain 

00 1 2 

(5-6) V 1 =-. 

This is equivalent to Euler's result £(2) = 7r 2 /6 because letting 

oo oo 1 

(5-7) c(*) = £~ and (o dd (k) = j2^— T] 

j= i j j= i y j > 

we have the following lemma. 
Lemma 5.2. Fork > 1, 

(5-8) ak) = 7^—;Co dd {k). 



2 k - 1 



Proof. 



n=l n odd n even 



oo ^ 
n=l v 7 

=C«h(*) + 2*C(*)- 

Now solve for C(^)- D 
Therefore by (15.61) . 

(5-9) C(2) = ~( odd (2) = y- 

Note that substituting z = e(t/N), < t < 1, into (I5.2p and applying the same manipu- 



lations as above gives 



TV 2 



n + t) 2 sin 2 (7rf)' 
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From this one can produce results analogous to Euler's, but where the sum is over an 
arithmetic progresion. 

In the next section we generalize the method to evaluate the zeta-function at the positive 
even integers. 



6. Zeta at the Even Integers 

The same calculations as above with a little more book keeping allow us to compute ((2n) 
for all positive integers n. This general result was known to Euler and there are many proofs 
of it. To state the result we need to introduce the Bernoulli numbers. The Bernoulli numbers 
Bk for k > 1 are defined as the coefficients of the Taylor series for 

z B k h 1 1 z 2 1 z 3 

6.1 => -^-z k = 1- -z + — - + ■••. 

v ; e z - 1 ^ k\ 2 6 2! 30 3! 

k>l 

For more about the Bernoulli numbers, see [2]. We can now state the general theorem 
which gives the values for ((2n). 

Theorem 6.1. For n > 1, 

1 (-l) n+1 2 2n - 1 n 2n B 2n 



(6.2) C(2n) := 



k>l 



k 2n " (2n) 



We provide a proof of Theorem l6.1l by differentiating the polynomials z N — 1 in two different 
ways. In each case we will have an expression involving recursively defined constants that 
are well known combinatorial objects. Using a couple of identities for these numbers and 
letting N tend to infinity we will be able to evaluate C(2^)- 

6.1. Repeatedly "Differentiating". As in the calculation of C(2), we will see that taking 
the logarithm of our polynomial before differentiating will allow us to express things in a 
useful way. However, unlike that calculation we will not simply take the derivative. Instead 
of repeatedly applying the operator j- we apply the operator z-^-. We will refer to this as 
"differentiation," where the quotes indicate that we are applying z-^ instead of ^. 
We wish to apply z4- many times to the function 



(6.3) log(^-l)= ^g(z-e 

j mod N 



To understand high "derivatives" of the right hand side we will consider the high "derivatives" 
of the function log(,2 — a), where a is a complex number. We have: 

Proposition 6.2. Define bk^ m by b\ t i = 1, and &o,m — for all m > 0, and 

(6.4) &fc, m +i = kb k , m - (k - l)b k -i, m 

for m > 1 and < k < m. We have 

(d\ m m 
z—j \og(z-a) = J2 b k, m z k (z-ay k . 
Z J k=i 
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Proof. Begin by noticing that 

(6.6) (z-^j {z k {z - a)- k ) = k (z k (z - a)~ fc - z k+ \z - a)^ 1 ) 

for k > 1. We will now establish the result by induction. It is true for m = 1. Assuming it 
is true for for all r < m, with m > 2, we have 

E lo g(^-a) = (^)X)6fc,m-i«*(«-ar* 

' j mod AT ^ ' k=l 

m—1 



b ^m-i {kz k (z - a)- k - kz k+ \z - a)-^) 

k=l 
m 

(kh,m-i ~(k- l)6jfe_i iro _i) z k (z - a)~ k 



k=l 
m 



^6 fejm 2; fc (z-a) fc . 



fc=i 



The first step follows by the induction hypothesis and the final step follows by the definition 
of the 6fe >m . □ 

Using this result we see that 

z Tz) E ^g{z-e{ 3 /N)) = J2b k , m z\z-e(j/N))- k 

j mod N j mod N k=l 



-k 



(6-7) = X>,m *{z-*U/N)Y 

k=l j mod N 

which is an expression for the "derivatives" of the right side of (16. 3p . Next we will find an 
expression for the derivatives of the left side. 

Proposition 6.3. Define c k ^ m by = 1, c^i = for all k > 1, and c 0)m = for all m > 0, 

and 

(6.8) c Km+ i = - (kc kyTn + (m + 1 - fc)c fc _ lim ) 

/or > 1 and m > 1 . VFe /iai>e 

Proof. Define Q m ( z ) '■= (z N ~ l) m ( z ^) m l°S — l) • An easy inductive argument shows 
that Q m {z) is a polynomial. In the notation of the Proposition we wish to prove that 

m 

(6.10) Q m {z)=N m Y,Ck, m z Nk . 

k=l 

We will obtain a recursion for Q m and then using that recursion we will find a recursion for 
its coefficients, c km . 
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By definition we have 



Q m +i(z) _( d\ Qm(z) 



(z N - l) m+1 \ dzj (z N - l) m 

_ z(z N - l)Q' m (z) - mNz N Q m (z) 
{z N - l) m+1 

Thus we obtain the recursion 

Q m+1 (z) =z(z N - l)Q' m (z) - mNz N Q m (z) 

=z N+1 Q' m (z) - zQ' m (z) - mNz N Q m (z). 

By induction we have the following 

m 

z N+1 Q'(z) =N m+1 Tkc k . m z N ^ 



k=l 
m 

zQ' m (z) =N m+1 J2kc k 



k=l 

Nk 



z 

k=l 



mN Z N Q rn {z) =N m+1 J2™Ck,mZ N{k+1) 

k=l 

Combining these we obtain the recursive formula for the constants, Ck, m - D 

6.2. Combinatorial Numbers. In Propositions 16.21 and 16.31 we introduced two sets of 
numbers {&fc,m} and {cfc jm }, and those numbers satisfied particular recursion formulas. It 
turns out that those numbers are related to other famous numbers from combinatorics: 
Stirling numbers of the second kind and the Eulerian numbers. 

The Stirling numbers of the second kind S(n, k) count the number of ways to partition a 
set of n elements into k nonempty subsets. They can be described recursively by S(l, 1) = 1 
and the recursion 

(6.11) S(m, k) = S(m - 1, k - 1) + kS(m - 1, k). 

On the other hand, the Eulerian number efc, m is the number of permutations of the numbers 
1 to m in which exactly k + 1 elements are greater than the previous element. The Eulerian 
numbers are given recursively by ei t \ = 1 and e^i = for k > 1, and 

(6.12) e kym = {m-k + l)e fc _i ;m _i + ke k ^ m _ x . 

It is the recursive definitions which we will make use of. 
We have the following: 

Proposition 6.4. In the notation above, 

b k , m =(m - l)\S(m, k) 

Ck^m ( 1) &k,m- 

In particular, b mm = (— l) m_1 (m — 1)!. 
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Proof. We will not prove the first of these two results because the only fact we need about 
b nm is b m ^ m = (— l) m-1 (m — 1)!. To prove this, it is enough to note that by the recursion in 
Proposition 16.21 we have 



(6.13) b mtm = TO& mjm _l - (iTL - l)6 m _l, m _l = -(m ~ 1)6. 



m— l,m—l ■ 



Thus the result follows by induction and the fact that b\ t i = 1. 

By the definition of Ck, m in Proposition 16. 31 we see that Ck, m and (— l) m_1 efc im _i satisfy the 
same recurrence relation with the same initial conditions. □ 



The following relates the Eulerian numbers that appear in our formulas to the Bernoulli 
numbers which appear in the formula for C(2n). This Lemma is Exercise 72 in Chapter 6 
of0. 

Lemma 6.5. 



(6.14) £(-l) fc c M = (-l) n 2 n (2 n - 1)^. 

fc=i 

6.3. Limit as iV — > oo. We now return to our discussion of the two expressions for the 
"derivatives" of \og(z N — 1). As in the calculation of C(2), to evaluate C(2n) we begin by 

equating the two expressions for (z-^)^ \og(z N — 1) given by equation (16.71) and Proposi- 
tion [631 Using the Taylor series expansion that we used in the previous section and similar 
analysis we are led to the evaluation of Codd(2n). 



Proof of Theorem \6.1[ Set m = 2n. Beginning with equation (16.71) and setting z = e(l/2N) 
we have 

(7 \ 771 777 
z te) E ^g(z-e(j/N))\ z=e{1/2N) =J2b k , m Yl e{l/2N)\e{l/2N)-e{j/N))- k 
j mod JV k=l j mod TV 

(2j-i; 



k=l j mod N 



e 



2N 



As in the calculation from Section \5\ we will use the Taylor series 



(6.15) (l_e(a;))- 1 = -^ + 0(l), 

2mx 
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which is valid for < \x\ < 1. To use the Taylor series we must have —N/2 < j < N/2 and 
so we obtain 



k=l -N/2<j<N/2 



2N 



—k 



k=l -N/2<j<N/2 



2N 



b V 

J m,m / j 



-N/2<j<N/2 



2m(2j - 1) 
2N 



O 



N 



k-l 



(2j - l)*-> 
2H23-T) ) + -( JV "- 1 ) 



m— 1 



fc=2 -N/2<j<N/2 



2N 



(6.16) 



+ E 

-N/2<j<N/2 

7ii 



2N 



2m(2j - 1 
(2j 



27Ti(2j - 1) 

+ 0(1 



+ O k 



N 



k-l 



(2j - !)*-> 



E 

N/2<j<N/2 



^y) + O rn (N™- 1 ) + 0(AMog(AT)). 



We write O m because the implied constant depends on m. The O(NlogN) term is only 
necessary when m = 2. 

Next we will turn to our other expression for (z^) log(z N — 1), which was given in Propo- 
sition [6731 As before, evaluating at e(l/2N) we have 



^ m Er=i^ m e(l/2AQ 
(e(l/2N) N - l) m 



Nk jym "' 



ZXm(-i)' 



(6.17) 



v ' k=l 

N m ( B 
' _]^ m 2 m (2 m 1) 

m 



where the second equality follows from the Lemma 16.51 and we use that m is even. By 
equation (I6.3P we have that equations (16. 16[) and (16.171) are equal. Setting the last lines of 
each equation equal and dividing by N m we have 
(6.18) 



(2 m -l)— ) =b m . m (mY 
m 



£ ((27^1)) + ^{N- 1 )+0(N 1 - m \og(N)). 



-N/2<j<N/2 



Letting N go to infinity we obtain 



(6.19) 



oo ^ 



(2 m - l)B ri 



(«) m J f^ 00 (2j-l) r 



m 
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To finish the proof we use the facts that b m ^ m = (— l)' m_1 (m — 1)! from Lemma 16.21 and 
Lemma [5.21 to obtain 



C(m) 



l-2 m ) ml 
2 m 1 ] 



Using m = 2n we have 



(1 - 2 m ) 2 .f-^ (2j - l) m 

2 m-i (2 m - l)( 7 ri) m J B m 
~(1 - 2 m ) ml ' 

(6.20) C(2n) = 1 J (2n) , ^, 

the desired result. □ 

Instead of the substitution z = e(l/2N) one can set z = e(A/BN) to obtain an expression 
for 



(A + Bnf 

This can then be used to evaluate Dirichlet L-functions L(s,x), with character x, at the 
positive even (odd) integers if x is even (odd). 
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